Broken symmetry and the variation of critical properties in the phase behaviour of 

supramolecular rhombus tilings 
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The degree of randomness, or partial order, present in two-dimensional supramolecular arrays 
of isophthalate tetracarboxylic acids is shown to vary due to subtle chemical changes such as the 
choice of solvent or small differences in molecular dimensions. This variation may be quantified 
using an order parameter and reveals a novel phase behaviour including random tiling with varying 
critical properties as well as ordered phases dominated by either parallel or non-parallel alignment 
of neighbouring molecules, consistent with long-standing theoretical studies. The balance between 
order and randomness is driven by small differences in the intermolecular interaction energies, which 
we show, using numerical simulations, can be related to the measured order parameter. Significant 
variations occur even when the energy difference is much less than the thermal energy highlighting 
the delicate balance between entropic and energetic effects in complex self-assembly processes. 



Two-dimensional molecular networks provide an at- 
tractive and highly flexible route to the formation 
of surfaces with specific structural arrangements and 
functionalities-^. Within the broad range of structures 
that have been explored, random molecular networks 
have gained considerable attention recently as exemplars 
of two-dimensional, nonperiodic glassy systems^—. Al- 
though these systems are of great interest, and represent 
examples of complex self-assembled sytems, a quantita- 
tive measure of the degree of randomness present in a 
molecular array, which is required for a theoretical un- 
derstanding of their properties, is not generally available. 
In this paper we report a set of supramolecular arrays 
in which the randomness varies due to small changes in 
chemical environment or molecular geometry. We define 
an order parameter which quantifies these variations and 
provides a systematic comparison between supramolecu- 
lar arrays prepared using similar methodologies. Further- 
more, we show that the observed variations may be ratio- 
nalised within a complex phase diagram which can be re- 
lated to the adsorption of dimers on a periodic lattice 7 ^, 
a paradigmatical problem in statistical mechanics which 
has been studied for several decades, and is a well-known 
example of a physical system which can be mapped onto 
a non-periodic tiling of the planed— . In particular, us- 
ing numerical simulations, we relate the order parameter 
to an orientational dependence of the intermolecular in- 
teractions which stabilise the network. 

We focus on molecular rhombus tilings formed at 
the liquid-solid interface following the deposition of sol- 
vated molecules of interest on highly-oriented pyrolytic 
graphite (HOPG) substrates. In recent work we have 
shown^ that such arrays support the formation of a ran- 
dom tiling, the first experimental realisation of the dimer 
adsorption problem. The molecular networks formed 
on graphite are imaged using scanning tunnelling mi- 
croscopy (STM) and the observed structures can be di- 



rectly mapped onto rhombus tilings. In-plane stabil- 
isation of the molecular networks arises from hydro- 
gen bonds formed between carboxylic acid groups ^ 17 ' 18 . 
For the tetracarboxylic acids considered here this re- 
sults in two orientations of intermolecular bonding with 
the molecular backbone axes either parallel or non- 
parallel at 60°, as shown in Fig. la. To form a 
random tiling a close match between the two dimen- 
sions marked d\ and di in Fig. la is required. Ex- 
amples of suitable molecules are p-terphenyl-3,5,3",5"- 
tetracarboxylic acid (TPTC) (Fig. lj), with d\ = 8.7 A 
which we have studied previously^, and 1,4-diphenyl- 
l,3-butadiyne-3,3",5,5"-tetracarboxylic acid (DPBDTC) 
(Fig. li), with di = 9.5 A synthesised specifically for this 
study; d 2 = 9.6 A for both molecules. Several closely re- 
lated molecules have been studied previouslj*^ 7 - - — , but 
do not display an extended random tiling, although frus- 
trated crystallisation has been observed^. To tile the 
plane perfectly a rhombus must have internal angles of 
60° and 120°, a lozenge^, for which di = d^. 

A random tiling is considered ideal when the energies 
of the non-parallel and parallel alignments, given by En 
and ep respectively, are degenerate, En = ep^—. In 
this case the tiling configuration is determined solely by 
maximising configurational entropy. Such tilings possess 
no translational order and display correlations in tile ori- 
entation which have a logarithmic dependence on spatial 
separation, a signature of a 'Coulomb' phased. However, 
if there is an energetic asymmetry, A = e^v — sp =/= 0, the 
degeneracy between the exponentially large number of 
tilings is broken. In this case an additional internal en- 
ergy contribution to the free energy competes with the 
entropic term (rhombus tiles with asymmetric interac- 
tions are described as interacting in many theoretical 
models^—). For small A, random tiling phases (with 
properties which vary continuously with A) are expected, 
but large A values are predicted to lead to transitions to 
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FIG. 1: Tetracarboxylic acid supramolecular assemblies and rhombus tilings, a, Parallel and non-parallel intermolec- 
ular bonding orientations, with backbone dx, and bond, cfe, lengths indicated, overlaid onto the corresponding rhombus tile 
representations, b, Schematic phase diagram of the interacting rhombus tiling model as predicted by theory. Three phases 
are expected depending on the interaction energy A: a random tiled phase with critical correlations for small |A| (a so-called 
"Coulomb" phase), an ordered phase dominated by non-parallel bonding at large positive A, and a second ordered phase dom- 
inated by parallel bonding at large negative A. The transition between the random tiled and non-parallel phases is expected 
to be continuous at T > (indicated by a full line in the sketch), and that to the parallel phase first-order (dashed line), 
c-h, STM images (sections of larger scans) of tetracarboxylic acid supramolecular networks at alkanoic-acid-HOPG interfaces 
and the corresponding rhombus tilings, in order of decreasing ip. Experiments are labelled I-IV and are performed at room 
temperature (apart from IV which was at 60°C) using different combinations of TPTC, DPBDTC, coronene and solvents as 
specified in text boxes. STM image contrast originates from molecular backbones (and, in h, coronene). All scale bars 50 A, 
see Methods for imaging parameters, i-k, Molecular ball and stick diagrams of; i, DPBDTC; j, TPTC; and, k, coronene. 1, 
Diagram of coronene adsorbed at the vertex of six TPTC molecules. 



ordered phases^—. The parameter A may be consid- 
ered to be analogous to the coupling constant J which ap- 
pears in Ising models of magnetism and leads to ordered 
phases with broken symmetry for | J| > fcgT, which may 
be either ferromagnetic or antiferromagnetic depending 
on the sign of J, while entropic terms in the free en- 
ergy dominate for \J\ << fc^X 1 in which case a param- 
agnetic phase occurs. Our experiments and simulations 
show that this rich phase behaviour may be investigated 
by preparing molecular arrays with differing values of A. 

Figures lc to lh show images of molecular net- 
works prepared under several different conditions (see 



Methods), referred to as Experiments I to VI, respec- 
tively, and the corresponding rhombus tiling representa- 
tions, where each molecule is represented by a rhombus 
coloured according to its orientation. To characterise 
our experimental tilings we define an order parameter, 
ip — (nop — pon) / (nop + pan) , where n and p represent 
the fraction of rhombus tile junctions in non-parallel and 
parallel orientations, respectively, and hq ~ 0.608 and 
Po ~ 0.392 are the equivalent values for a defect-free, 
ideal, random tiling and are estimated numerically (see 
Methods). As such, i/j = 1 in a fully parallel phase, 
ip = — 1 in a fully non-parallel phase, and ip = for an 
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ideal tiling, ip is calculated for each Experiment (I to 
VI), and the images (Figs, lc to lh) are placed in order 
of decreasing ip, spanning ip — 0.22 to ip = —0.68. 

Fig. lb is a schematic of the expected equilibrium 
phase diagram-^— . There are two relevant thermody- 
namic parameters: temperature, T, and energetic bias, 
A. Ideal random tilings are observed for A = 0. For 
A > non-parallel bonding is favoured, resulting in in- 
creasingly negative values of ip. For sufficiently large A 
the system orders into a crystalline phase dominated by 
non-parallel bonds. This transition occurs at A/fcgT = 
0.454(3)^£. For A < parallel bonding is favoured, 
resulting in increasingly positive values of ip. The ran- 
dom tiled phase extends until A/fc^T « —0.3 where 
the system undergoes a first-order transition^ to a crys- 
talline phase dominated by non-parallel bonds. We re- 
iterate the analogy with magnetic systems where tran- 
sitions from disordered to ordered phases occur as the 
absolute value of the coupling constant is increased be- 
yond some critical value. 

It is clear from Fig. 1 that different experimental con- 
ditions give rise to tilings with varying degrees of order. 
In Experiment I we investigated DPBDTC in nonanoic 
acid (Fig. lc). For this arrangement there are more par- 
allel bonds than expected for an ideal random tiling and 
ip = 0.22. A solution of TPTC in heptanoic acid (Exp. 
II, Fig. Id) results in a molecular tiling with ip — —0.08. 
This is the molecular tiling closest to ideal (i.e. ip = 0) 
of all the systems we have studied. For TPTC a change 
of solven t 21 i 22 from heptanoic acid to octanoic (Exp. Ill, 
Fig. le, ip = —0.25) or to nonanoic (Exp. V, Fig. lg, 
ip = —0.43) acid, solvent molecules with similar chem- 
istry but slightly larger size, leads to a progressive in- 
crease in the fraction of non-parallel bonds. Experiment 
IV is performed with TPTC and nonanoic acid at an el- 
evated temperature of 60 °C (Fig. If, ip = —0.35). Com- 
paring this molecular tiling with the equivalent tiling pre- 
pared at 19 °C (room temperature) we see that the in- 
crease in temperature leads, as expected, to a value of ip 
closer to zero. Finally, the most negative order param- 
eter, ip = —0.68, is observed when coronene (Fig. Ik) is 
added to a solution of TPTC and nonanoic acid (Exp. 
VI, Fig. lh). Coronene becomes incorporated in the 
molecular network in a site which enhances the stabil- 
isation of the non-parallel arrangement (Fig. 11), a cap- 
ture process consistent with previous studies of coronene- 
carboxylic acid and other system o 19 i 23 ~ — . 

These results show that small modifications in experi- 
mental parameters lead to the exploration of the rhombus 
tiling model phase space. To determine the placement 
of the experimental structures on the phase diagram in 
Fig. lb, we have undertaken extensive numerical simu- 
lations of rhombus tilings. The details of the numerical 
scheme are provided in the Methods section and repre- 
sent a generalisation of previous wor k 26 i 27 to the case 
A ^ 0. Fig. 2a shows the numerically computed order 
parameter ip as a function of scaled bias A/ksT. Verti- 
cal dashed and solid lines indicate where transitions from 



random to ordered phases occur. The six experiments are 
marked on the numerical curve according to their mea- 
sured ip values. Experiments II to V span the random 
tiled region with negative ip; Experiment VI lies beyond 
the phase boundary in the non-parallel phase; Experi- 
ment I is close to the phase boundary with the ordered 
parallel phase. 

The numerical results in Fig. 2a provide a relationship 
between the order parameter ip and A/fc^T which allow 
an estimation of the latter for each of our experiments. 
Furthermore, using the known temperature at which the 
experiments were performed we can estimate values of 
A, as tabulated in Fig. 2c. A comparison of Experiments 
IV and V is of particular interest: both were performed 
using TPTC in nonanoic acid but with molecular net- 
works prepared at different temperatures. Accordingly 
we would anticipate the same value of A and we find ex- 
cellent agreement between the values inferred from simu- 
lations. This result provides strong support for our anal- 
ysis and confirmation that the molecular arrays can be 
understood in terms of equilibrium structures. From our 
experiments it is clear that measurable variations in the 
tiling statistics may be identified even for variations of A 
of < 5 meV, much less than the thermal energy at room 
temperature. 

A second experimental measure, which can also be 
compared with numerical results, is the concentration of 
defects. These 'half-rhombus' vacancies are topological 
defects and are present in several of the images shown 
in Fig. 1. For an ideal tiling their concentration, c, has 
been shown to be temperature dependent^ (we choose 
an effective temperature in our simulations to match, 
approximately, the maximum defect density observed in 
experiments) and we show in Fig. 2b that c is also con- 
trolled by A. We highlight the non-monotonic depen- 
dence on bias which results from the strong suppression 
of isolated triangular defects in the ordered phases. The 
experimentally-determined defect data is plotted against 
the theoretical curve using the value of A/fc^T obtained 
from Fig. 2a. We note that the non-monotonic behaviour 
and suppression of defect formation in the ordered phases 
is reproduced in our data thus providing further support 
of the validity of our approach to determine the value of 
A for our experiments. 

Rhombus tilings are commonly analysed as a projec- 
tion of the interface of a simple cubic lattice onto a (111) 
plane 9 . In this approach each vertex is assigned an effec- 
tive height h{x, y) (x and y are the co-ordinates of the 
vertex) using a simple procedure (outlined in Fig. 3j). An 
analysis of our images using this approach shows explic- 
itly the difference between random and ordered phases as 
shown in Fig. 3. Here we show larger area images for Ex- 
periments I, IV and VI, in three different representations. 
The first is the tiling representation, as in Fig. 1 (Figs. 3a, 
3d, and 3g). The second representation is given by the 
height map h(x,y) (Figs. 3b, 3e, and 3h): Experiment I 
shows a clear height gradient, Fig. 3b, indicative of its 
closeness to the phase boundary to the parallel phase; 
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FIG. 2: Comparison of experiments to simulations of 
the rhombus tiling model, a, Order parameter ip as a 
function of scaled bias A/ksT. The full curve is the result 
of numerical simulations of the rhombus tiling model (see 
Methods) for ejv/fcgT = 5/3 (at this value the concentra- 
tion of defects c is comparable to that observed experimen- 
tally). The vertical lines indicate the approximate location of 
the continuous transition to the non-parallel ordered phase, 
A/k B T = 0.49(2) (slightly larger than that calculated^ for 
T = - see Methods and footnote-^), and the first-order tran- 
sition to the parallel phase, A/fcsT = —0.34(3). The experi- 
mental points are placed on the curve according to the mea- 
sured values of tp. h, Concentration of defects as a function 
of bias. The full curve is the numerical result. The abscissa 
values for experimentally measured defect concentrations are 
fixed by panel a: the experiments show the predicted non- 
monotonic dependence on bias, c, Estimate of bias A for 
the six experiments. Notice that A for Experiments IV and 
V, which correspond to the same system (TPTC in nonanoic 
acid) at different temperatures, is the same within experimen- 
tal error, as expected. 



in contrast the height maps for Experiments IV and VI, 
Figs. 3e and 3h, show averaged height gradients which 
are close to zero. This is expected for both the random 
and non-parallel ordered phase although these can be dis- 
tinguished by the greater height fluctuations present for 
Experiment IV consistent with a random tiling phase. 

The third representation illustrates the spatial vari- 
ation of non-parallel order. In a perfectly non-parallel 
configuration the unit cell (e.g. see schematic inset of 
Fig. lb) can be centred at any of the three sub-lattices 
of the triangular lattice. Choosing a sublattice is the 
symmetry breaking associated with the transition to the 
non-parallel phase and in general we would expect to ob- 
serve domains of non-parallel order corresponding to each 
of these three possible sub-lattices. Again a magnetic 
analogue is useful; in a ferromagnetic system domains of 
spin-up and spin-down magnetisation are formed whose 
size diverges as the system undergoes a transition from 
a paramagnetic to a ferromagnetic phase. For the case 
considered here there are three equivalent non-parallel or- 
dered states into which the system can crystallise, rather 
than the two states available for a simple magnetic sys- 
tem. In Figs. 3c, 3f, and 3i regions of non-parallel aligned 
tiles are coloured cyan/yellow/magenta depending on 
their assignment to one of the three possible non-parallel 
phases. For Experiment I, Fig. 3c, these domains are 
relatively small since the parallel tile alignment predom- 
inates. The domains are larger but finite for Experiment 
IV, Figs. 3f, in the random tiling phase. For Experi- 
ment VI, Fig. 3i, one domain has become dominant and 
is comparable with the size of the system, indicative of a 
transition to the non-parallel phase. 

Our results show conclusively that the phase space 
of two-dimensional rhombus tiles, which arises from 
asymmetric interactions, may be explored experimentally 
through the structural characterisation of molecular ar- 
rays. The changes in energy which give rise to different 
ordering are small compared to the thermal energy, an 
essential requirement for systems where entropy is signifi- 
cant. The dependence of the characteristic energetic bias, 
A, on solvent indicates that the effective energetic differ- 
ence between parallel and non-parallel molecular pairs 
adsorbed on graphite cannot be simply described by the 
interactions between two tetracarboxylic acid molecules 
and must include other more complex contributions. We 
highlight two possible influences of the solvent molecules 
on the effective interaction between adsorbed molecules. 
Firstly the co-adsorption of solvent molecules on the 
graphite surface is possible^, and expected, since the 
the network of tetracarboxylic acid molecules results in 
an array of nanopores which can accommodate guest 
molecules. As shown in a recent paper— these nanopores 
are hexagonal and may be classified into five inequivalent 
types which are formed from different combinations of 
parallel and non-parallel bonded molecular pairs. There 
is a significant variation of the capture of guest molecules 
in pores of different type^, and we therefore also ex- 
pect variations in the stabilisation energies of the co- 
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FIG. 3: Ordered phase analysis of tetracarboxylic acid supramolecular networks, a-c, Experiment I (DPBDTC 
in nonanoic acid): (a) tiling representation of a 50 nm STM image; (b) height field representation, showing that at this 
scale the system is 'tilted'; (c) non-parallel phase domain representation (we colour cyan/yellow/magenta the three possible 
different sub-lattice arrangements of the non-parallel ordered phase), d-f, Experiment IV (TPTC in nonanoic acid): (d) tiling 
representation (80 nm image), (e) height field representation; (f) non-parallel phase domain representation, e and f show that 
the structures of Experiment IV are neither tilted, nor display any significant non-parallel phase order, as expected of a random 
tiling, g-i, Experiment VI (TPTC and coronene): (g) tiling representation (100 nm image); (h) height field representation; 
(i) non-parallel ordered domain map. While h shows that the structure is overall flat, i shows that one ordered non-parallel 
domain has become prevalent indicating a spontaneous breaking of symmetry due to a transition to the ordered non-parallel 
phase, j, Scheme for calculating the height field: moving along the edges of tiles the field h(x,y) increases/decreases by the 
amounts indicated. The bar shows the scaling of height to contrast in the height representations; the same scaling is used in 
(b), (e) and (h). 



adsorption of solvent molecules in different pore types. 
Since the energy of co-adsorbed molecules makes a con- 
tribution to the overall stability^ of a given pore, this ef- 
fect provides a mechanism to break the degeneracy of the 
different pore types. Since different pores are formed by 
different combinations of parallel and non-paralell pairs, 
variations in guest inclusion stability will result in small 
effective energetic differences between the parallel and 
non-parallel molecular orientations. A second considera- 
tion is that the solvents provide an effective medium with 
varying dielectric constant, as pointed out for these par- 
ticular molecules by Kampschulte et.al^, which would 
give rise to a solvent dependence of any energetic im- 
balance mediated, for example, by van der Waals in- 
teractions, or hydrogen-bonding, between neighbouring 
molecules (see Cook et.ali^fi for a discussion of solvent 
effects on intermolecular interactions in a solution envi- 
ronment). 

A complete microscopic description of energy differ- 
ences on this scale are beyond the confidence levels 
normally associated with modern numerical approaches 
to molecular modelling (e.g. density functional theory, 
molecular dynamics). Furthermore the inclusion of sol- 
vent molecules makes a detailed microscopic description 
and prediction of energies particularly problematic. How- 
ever, our experiments give clear limits on the tolerance 
of the random tiling phase to symmetry breaking - a dif- 



ference in binding energies of > lOmeV will drive the 
structures into an ordered phase. 

Our study demonstrates that prototypical self- 
assembled systems such as supramolecular networks can 
reveal a richness of exotic phases which have been stud- 
ied widely by theorists and are analogous to those found 
in complex condensed-matter systems. In addition our 
work highlights the role of entropy in the balance be- 
tween order and disorder in molecular templates. 



Methods 

Chemicals. l,4-diphenyl-l,3-butadiyne-3,3" ,5,5"- 
tetracarboxylic acid (DPBDTC) and p-terphenyl- 
3,5,3" ,5"-tetracarboxylic acid (TPTC) were synthesized 
in-house, details of their synthesis can be found in the ac- 
companying Supplementary Material and Ref. H] respec- 
tively. Other chemicals used were coronene (>95.0%, 
Fluka), heptanoic acid (>97%, Sigma), octanoic acid 
(>98%, Aldrich), and nonanoic acid (>95%, Fluka). 

STM experiments. STM images were acquired with 
an Agilent Technologies 4500 PicoPlus STM using a Pi- 
coScan controller and STM tips formed from mechani- 
cally cut Ptlr (80:20) wire. Saturated solutions of the 
desired adsorbate (DPBDTC or TPTC) were prepared 
by placing an excess of solid in the desired solvent (hep- 
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tanoic, octanoic, or nonanoic acid). For Experiment VI, 
the solution used combined equal volumes of a satu- 
rated solution of TPTC in nonanoic acid and coronene 
in nonanoic acid (cone. 1.5 x 10~ 4 mg/mL). All so- 
lutions were ultrasonically agitated to ensure complete 
dissolution. To form the supramolecular networks, a 
10 fjlj droplet of solution was deposited onto a mechani- 
cally cleaved HOPG substrate. Imaging at the alkanoic- 
acid-HOPG interface commenced immediately after ap- 
proach of the STM tip. For Experiment IV, the substrate 
was heated to 60 °C prior to solution deposition and 
held at this temperature for the duration of the exper- 
iment. Figure 1 STM imaging parameters (Vu p /I): lc, 
+1.1 V/10 pA; Id, +1.2 V/15 pA; lc, +1.2 V/5 pA; If, 
+1.25 V/12 pA; lg, +1.0 V/13.5 pA; lh, +1.0 V/30 pA. 

Numerical simulations. The rhombus tilings were 
simulated using the same model as in Ref. [26| sup- 
plemented by interactions which varied depending on 
whether neighbouring tiles were parallel (binding energy 
ep) or non-parallel (binding energy En)- Simulations 
were performed with continuous time Monte Carlo which 
allows to access very long times at low temperatures (i.e. 
very low number of defects) on systems sizes between 
N = 18 2 to 99 2 and averaged over 10 3 samples for each 
state point. The location of the continuous transition 



from the random tiled to the ordered non-parallel phase 
at T ^ (the location at T = is known exactly^) was 
determined from finite size scaling analysis of susceptibil- 
ities and higher order cumulants of ip and of the 'colum- 
nar' order parameter (difference in density of non-parallel 
order in each sublattice) as in Ref. [ID; the location of the 
first-order transition to the ordered parallel phase from 
the jump in ip and the 'tilt' order parameter (difference 
in density of each kind of tile) , cf. Ref. [TBI 
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